Leading neutrons from polarized pp collisions 
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Abstract. We calculate the cross section and single-spin azimuthal asymmetry, A„{t) for inclusive 
neutron production in pp collisions at forward rapidities relative to the polarized proton. Absorptive 
corrections to the pion pole generate a relative phase between the spin-flip and non-flip amplitudes, 
which leads to an appreciable spin asymmetry. However, the asymmetry observed recently in the 
PHENIX experiment at RHIC at very small \t \ ~ 0.01 GeV^ cannot be explained by this mechanism. 

Keywords: neutrons, pions, polarization, absorption 
PACS: 13.85.Ni, ll.80.Gw, 12.40.Nn, ll.80.Cr 



INTRODUCTION 

We are considering the reaction of inclusive neutron production pp -^ Xn at large 
Feynman Xf as was measured recently by the PHENIX experiment at RHIC [1] at 
■\/s = 200 GeV. The preliminary data [2] are shown in Fig. 1. The measurements were 
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FIGURE 1. Single-spin asymmetry A„ in reaction pp -^ nX measured at y/s = 200GeV [2] 

done with a transversely polarized proton beam and the neutron was detected in the very 
near forward (or backward) direction relative to the polarized beam with a transverse 
momentum pj of the order of 100 MeV. An appreciable single spin asymmetry was 
found for forward neutrons, for events with neutron momenta larger than 30 GeV A^? = 

(— 9.0±0.6±0.9)% with an additional scale uncertainty f 1.0 _q 24 ) • Data also show 



no dramatic variation with xf- According to the Feynman scaling, one should not expect 
any strong energy dependence of Aa^(z). 

At the same time, neutrons produced with xf < demonstrate small asymmetry, 
consistent with zero. This fact is explained by the so called Abarbanel-Gross theorem [3] 
which predicts zero asymmetry for particles produced in the fragmentation region of an 
unpolarized beam. This theorem was proven within the approximation of triple-Regge 
poles. Regge cuts breakdown this statement, but the corrections calculated in [4] turned 
out to be rather small, less than 1%. 

The cross section of of neutron production in pp collisions was calculated recently 
in [5] where it was found that previous calculations [6, 7], although agreed with the 
ISR data, grossly underestimated the magnitude of absorptive corrections. New mea- 
surements done recently show that the normalization of the ISR data is overestimated 
by factor two. This was observed in pp collisions in the NA49 experiment [8], and in ep 
[9] collisions. Here we concentrate on the spin dependence of neutron production in the 
fragmentation region of the polarized beam caused by the absorptive corrections. 

ABSORPTIVE CORRECTIONS TO THE PION POLE 

Pion exchange contribution to the amplitude of neutron production, pp -^ nX in Bom 
approximation has the form 

VZ 

where d are Pauli matrices; i^p „ are the proton or neutron spinors; qj is the transverse 
component of the momentum transfer; q^ = {I —z)mf^. Fractional momentum z is 
related to the invariant mass Mx of X, z ^ 1 — M^/s for 1 — z <^ 1, where ^/s is 
the cm. energy of pp collision. The 4-momentum transfer squared t has the form, 

t = -{.ql+<iT)/z- 

The differential cross section of inclusive neutron production reads [11, 12, 5], 

da^ _ I .. .|2_ 

= (^)'k|G2..,„W|r?.(0P(l-^)^-''^'^«c7,^^(^' = M|), (2) 

where rinit) is the signature factor, and G^+ Jyt) = exp(?7?j) is the pionic formfactor 
(see details in [5]). The results of calculation with the Bom approximation, Eq. (2), at 
^/s = 200, 62.7 and 30.6GeV are depicted versus ISR data [13] in Fig. 2. The data 
are shown at two energies ^/s = 30.6 GeV and 62.7 GeV. Correspondingly, we used 
these energies in our calculations. One can see that the Born approximation considerably 
exceeds the data and leads to a wrong energy dependence. 

Absorptive corrections, or initial/final state interactions look simpler and factorize in 
impact parameters. Therefore, we switch to a Fourier transformed amplitude (1) which 
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FIGURE 2. Left panel: Born approximation (dashed curve) for neutron production and ISR data [13] at 
V^ = 60.7 GeV and pj — 0. Two solid curves upper and bottom show the effect of absorptive corrections 
calculated in the dipole approach (x5'^''0 and in hadronic representation (x5''""^'0 respectively. Right 
panel: Energy dependence of inclusive neutron production. Three upper curves present the forward cross 
section at ^s — 30.6GeV (solid), 62.7GeV (dashed) and 200GeV (dotted-dashed) calculated in the Born 
approximation. These curves corrected for absorption ( x S'^''>') are presented by three curves at the bottom. 
Data at ^s ~ 30.6GeV and 62.7 GeV [13] are depicted by squares and triangles respectively 



has the form, 
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Here 



0(f(^^) 






a3gi^e^{b,z)-i^ef{b,z) 



(fqe^^<^^{qT. 



^(^)^■^^o(Vi8) + Y3^ [Ueb)-Ub/m 



(3) 



(4) 



Of{b,z) = \ld\e'^Hb-q)^\qT,z) 



and 



TTiH 1 



eKi{eb)--^Kiib/P) 



N{z) = \gn^pnz{l-z)<^-"^^^i/'^e-'''^l/^A,,^x{M^x 



X) 



£^ = ql + zml, 
1 



r = -[i?f-<ln(l-z)]. 

z 



(5) 



(6) 



The process under consideration at large z — > 1 is associated with creation of a 
large rapidity gap (LRG), Ay = |ln(l — z)|, where no particle is produced. Absorptive 
corrections may be also interpreted as a suppression related to the survival probability 
of LRG. 

To identify the projectile Fock states responsible for absorptive corrections we start 
with Fig. 3a containing the amplitude of pion-proton inelastic collision K + p ^*X. The 
latter is usually described as color exchange leading to creation of two color octet states 
with a large rapidity interval ~ \n{M^/sQ) {sq = IGeV^) as is illustrated in Fig. 3b. 
Perturbatively, the interaction is mediated by gluonic exchanges. Nonperturbatively, e.g. 
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FIGURE 3. a: Born graph with single pion exchange and excitation of the projectile proton, p + n^X; 
b: inelastic proton-pion interaction interaction, p + n ^ X, via color exchange leading to production 
of two color-octet dipoles which hadronize further to X; c: Fock state representation of the previous 
mechanism. A color octet-octet dipole which is a 5-quark Fock component of the projectile proton, 
interacts with the target proton via 7r+ exchange. This 5-quark state may experience initial and final state 
interaction via vacuum quantum number (Pomeron) exchange with the nucleons (ladder-like strips). 

in the string model, hadron collision looks like crossing and flip of the strings. 

According to Fig. 3b one may think that this is the produced color octet-octet state 
which experiences final state interactions with the recoil neutron. On the other hand, 
at high energies multiple interactions become coherent, and both initial and final state 
interactions must be included. This leads to a specific space-time development of the 
process at high energies, namely, the incoming proton fluctuates into a 5-quark state 
|{3^}8{^^}8) long in advance of the interaction with the target via pion exchange, as is 
illustrated in Fig. 3c. 

This picture can be also interpreted in terms of the Reggeon calculus [5], which helps 
to identify what was overlooked in the calculation of absorptive corrections done in [6]. 

The resulting amplitude exposes both initial and final state attenuation of the 5-quark 
state. 



fp^nib.z)=f%,ib,z)Sib,z) 



(7) 



which are incorporated via the survival amplitude S{b,z). 

We evaluate the survival amplitude S{b,z) within two models based on the color- 
dipole and hadronic representations. 

Color-dipole model. 
According to the dual parton model [14] one can present the survival amplitude of a 
5-quark state with an accuracy I/N]: as. 



5(5?) (&) = S^'^'i\b)S^'i^\b) =\\- imr(3«)^(Z7)l [1 - Imr(«^)^(Z7) 



(8) 



similar to Eq. (11), The elastic amplitudes r^^^>P{b) of a color {33} dipole interacting 
with a proton can be calculated in terms of the partial dipole elastic amplitudes, for 
which a saturated model was proposed recently [15, 5]. The results for S^^'^' (b), S^'^^' (b) 
and their product S^^^^ (b) at ^ = 60GeV and z = 0.8 are depicted in Fig. 4 (left panel). 



0.75 - 



0.5 



0.25 



- Vs=60 GeV 


^^^^ 


z=0.8 


-'-^'^/'''^'"'^ 




■ ''^ y^ 




^/ X 


y 




/' 


/ 


/' 






/ 


/ / / 


^^\ s'^"' 


y' , / 








/ / 


^ s*'"' 


/ / 




/ / 




/ / 




/ / 




1 


1 1 



0.75 



0.5 



0.25 



- Vs=60 GeV 








z=0.8 


'''^ ''^^^^^^^ 




y y^ 


/* \ 




,-■ 


\/ 


/' / 


Y 


/' ^\/ 


/ \ c(^P} 


/" ^ A 




/" / / 




/* / / 


\ s'pp) 


/* / / 




/' ^ / 




/ / / 




— "'' ^ / 




y / 




_ - ^ / 




1 1 


1 



0.5 1 1.5 

b(fm) 



0.5 1 1.5 

b(fm) 



FIGURE 4. Partial survival amplitude S{b,z) at y/s ~ 60GeV and z — 0.8. Left panel: Survival am- 
plitudes 5P^)(fe,z) and S^^'^){b,z) depicted by dot-dashed and dashed curves, respectively. Solid curve 
corresponds to S^^^){b^z)- Right panel: Survival amplitudes evaluated in hadronic representation. 



Hadronic representation. 
One can expand the 5-quark Fock state over the hadronic basis, 

|{3^}8{^^}8) = do\p) +di\N7t)+d2\N27t) + ... . 



(9) 



Since the admixture of sea quarks in the proton is small, the projection of the 5-quark 
state to the proton must be small too. It is natural to assume that the amplitude di is the 
dominant one, since both states \{3q}s{^(l}s) and |A'';r) have the same valence quark 
content. Correspondingly, the absorption corrected partial amplitude gets the form. 



fp^n{b,z)=f'{b.z)S^^"^^\b). 



(10) 



where 



S^^'^''\b)=S''P{b)SPP{b) = [1 - lmYPP{b)] [1 - lmY''P{b)] . (11) 

The elastic partial amplitudes lraT^P{b) were extracted in a model independent way 
directly from data on elastic scattering [16, 5]. 

Our results for S^^^^^^' {b, z) are depicted in Fig. 4, right panel. The survival amplitude 
suppression factor S^^"''^''\B,z) is rather similar to one calculated in the dipole model, 

S^^'\b,z). 



ABSORPTION CORRECTED SPIN AMPLITUDES 



An example of absorption corrected impact-parameter-dependent spin amplitudes is 
shown in Fig. 5 for y^y = 60GeV and z = 0.8 in comparison with the Born approxi- 
mation. 
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FIGURE 5. Partial spin amplitudes for neutron production. Solid curves show the result of Born 
approximation. Dashed and dot-dashed curves include absorptive corrections calculated in the dipole 
approach (xS^^'<''{b,z)) and in hadronic model {xS^^ \b,z)), respectively. 

Now, it is straightforward to Fourier transform the absorption corrected amplitudes 
(7) back to momentum representation, 
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The forward neutron production cross section corrected for absorption is compared 
with the ISR data [13] in Fig. 2. The two models for absorption, dipole and hadronic. 



are presented on the left panel by solid curves, upper and bottom respectively. The 
results of both models are pretty close to each other, but substantially underestimate 
the ISR data. It was argued, however, in [5] that the normalization of these data is 
twice overestimated. This is supported by recent direct measurements of pp — > nX by 
the NA49 experiment [8] and by comparison with neutron production in DIS [9] (see 
discussion in [5]). Otherwise, the shape of the z-dependence is reproduced pretty well. 
Also the energy dependence of the cross section is much improved after inclusion of 
absorption. Apparently, the steep rise of the cross section with energy observed in Bom 
approximation is nearly compensated by the falling energy dependence of LRG survival 
amplitudes. The results for qt dependent differential cross section can be found in [5]. 
As far as the spin amplitudes are known, one can calculate the single spin asymmetry. 
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Our predictions for A]\/{qT,z) are presented in Fig. 6. Notice that the spin asymmetry 




FIGURE 6. The curves in the first and second left panels correspond to the absorption corrections cal- 
culated in the dipole and hadronic models respectively. The spin asymmetry is calculated at five different 
fixed angles between the neutron and initial proton momenta in the cm. of collision Oc.m. — 0.001 — 0.005. 
Correspondingly, neutron transverse and longitudinal momenta correlate, qj = 9c.m. z V*/2. The last panel 
shows the qr dependence of A^^ calculated in the dipole model at four fixed values of z = 0.6 — 0.9. 

depicted in the two left panels, demonstrate practically no model dependence of the spin 
asymmetry. 

Unfortunately, our calculations at transverse momenta as small as was measured at 
RHIC [1] considerably underestimate the data. This is not a surprise, q^ ~ O.OlGeV^ 
is so small, that hadronic spin effects in other reactions are usually vanishingly small. 
The spin asymmetry which arises due to pion pole with absorptive corrections, develops 
a double zero at qj -^ 0. This is because both, the spin-flip amplitude and the relative 
phase, Eq. (16), vanish. However, at larger transverse momenta the single spin asym- 
metry reaches an appreciable value. One should include other mechanisms, for instance 



interference of pion and ai exchanges, which is enhanced by a proper phase shift, n/l 
and energy independence of the diffractive amplitude Kp ^ a\p. 
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